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ABSTRACT 

Recent observations have shown that circumstellar and circumbinary discs in young 
stellar binaries are often misaligned with respect to the binary orbital plane. We ana¬ 
lyze the tidal truncation of such misaligned discs due to torques applied to the disc at 
the Lindblad resonances from the tidal forcings of the binary. We consider eccentric 
binaries with arbitrary binary-disc inclination angles. We determine the dependence 
of the tidal forcing strengths on the binary parameters and show that they are compli¬ 
cated non-monotonic functions of eccentricity and inclination. We adopt a truncation 
criterion determined by the balance between resonant torque and viscous torque, and 
use it to calculate the outer radii of circumstellar discs and the inner radii of circumbi¬ 
nary discs. Misaligned circumstellar discs have systematically larger outer radii than 
aligned discs, and are likely to fill their Roche lobes if inclined by more than 45° —90°, 
depending on the binary mass ratio and disc viscosity parameter. Misaligned circumbi¬ 
nary discs generally have smaller inner radii than aligned discs, but the details depend 
sensitively on the binary and disc parameters. 
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1 INTRODUCTION 

Protoplanetary discs are often found in binary systems. In 
general, because of the complex (and often turbulent) for¬ 
mation processes of stars, binaries, and discs (e.g.. Bate, 
Bonnell & Bromm 2003; McKee & Ostriker 2007; Klessen 
2011; Fielding et al. 2015), both circumstellar and circumbi¬ 
nary discs are likely formed with inclined orientations with 
respect to their binary orbital planes. Indeed, a number of 
binary young stellar objects (YSOs) are observed to con¬ 
tain circumstellar discs that are misaligned with the binary 
plane (e.g., Stapelfeldt et al. 1998, 2003; Neuhauser et al. 
2009). The orientation of jets in several unresolved YSOs or 
pre-main-sequence binaries also suggest misalignment (e.g., 
Davis, Mundt & Eisloffel 1994; Roccatagliata et al. 2011). 
Recently, Jensen & Akeson (2014) showed evidence that at 
least one of the circumstellar discs in the young binary sys¬ 
tem HK Tan is inclined with respect to the binary plane 
by at least 30°. Imaging of circumbinary debris discs indi¬ 
cates that while some are aligned with the binary plane, 
others can have significant misalignment (such as 99 Her- 
culis, which is misaligned by at least 30°; Kennedy et al. 
2012a, 2012b). The pre-main-sequence binary KH 15D is 
surrounded by a processing circumbinary disc inclined with 
respect to the binary plane by 10 — 20° (Chiang & Murray- 
Clay 2004; Winn et al. 2004; Capelo et al. 2012), and in the 
FS Tau system, circumstellar discs appear to be misaligned 


with a circumbinary disc (Hioki et al. 2011). Although cir¬ 
cumbinary discs may align with the binary plane on rela¬ 
tively short timescales, circumstellar discs are expected to 
maintain misalignment over timescales comparable to their 
lifetimes (Foucart & Lai 2014). 

In recent years, many exoplanets have been found in bi¬ 
naries, including both S-type planets (orbiting a single mem¬ 
ber of the binary) using radial velocity measurements (e.g., 
Hatzes et al. 2003; Chauvin et al. 2011; Dumusque et al. 
2012), and P-type (circumbinary) planets using the transit 
method (e.g., Doyle et al. 2011; Orosz et al. 2012; Welsh et 
al. 2012, 2014; Kostov et al. 2014). While the orbits of most 
of the circumbinary planets are consistent with alignment 
with the binary orbit, Kepler-413b appears to be slightly 
misaligned (by about 2.5°), possibly as a result of forming 
in a misaligned circumbinary disc (Kostov et al. 2014). 

The sizes of the protoplanetary discs (the outer radius 
of a circumstellar disc and the inner radius of a circumbi¬ 
nary disc) in which planets in binaries form are affected by 
gravitational interactions between the disc and binary. The 
size of the disc is of interest because it places restrictions on 
the process of planet formation. For example, the process of 
giant planet formation may be jeopardized in circumstellar 
discs with significantly truncated outer radii (Jang-Condell 
2015). 

The size of a pressureless circumstellar disc is deter¬ 
mined by the “static”, non-resonant tidal force from the bi- 
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nary companion. In particular, the disc size is determined 
by the location at which orbits around the primary star be¬ 
gin to intersect one another, in the context of the circular 
restricted three-body problem (Paczyhski 1977; Paczyhski 
& Rudak 1980). In this framework, the outer edge of a cir- 
cumstellar disc cannot extend past the 2:1 commensurability 
(i.e., the period ratio of the binary and a test mass in the 
disc must exceed two), except when the binary mass ratio is 
extreme. Pichardo, Sparke & Aguilar (2005) extended this 
analysis to eccentric coplanar binaries, and found that the 
disc must always be smaller than 73% of the average Roche 
lobe radius. However, Papaloizou & Pringle (1977) demon¬ 
strated that even a very small viscosity can prevent orbit 
crossings and allow the disc to be larger than the radius 
determined by non-intersecting orbits. 

In addition to static tidal interactions, gas discs (both 
circumstellar and circumbinary) in binaries are also subject 
to resonant tidal forcings. The motion of the binary excites 
spiral density waves at Lindblad resonances, where the nat¬ 
ural epicyclic frequency k of the disc is commensurate with 
the binary orbital frequency Hb. As a result, a torque is ex¬ 
erted on the disc at each resonance (Goldreich & Tremaine 
1979, 1980). Artymowicz & Lubow (1994; hereafter AL94) 
applied the formalism of Goldreich & Tremaine to determine 
the sizes of circumstellar and circumbinary discs in coplanar 
binaries. In their approach, the gravitational potential of an 
eccentric binary is decomposed into many Fourier compo¬ 
nents, each applying a torque on the disc at Lindblad reso¬ 
nances. The possibility of gap opening and disc truncation 
by each potential component is determined by the balance 
between the Lindblad torque and the viscous torque. The 
main goal of this paper is to generalize the results of AL94 
to circumstellar and circumbinary discs in binaries with ar¬ 
bitrary mutual inclinations. 

The behavior of Lindblad torques in misaligned discs 
has been explored in several recent works. Lubow, Martin 
& Nixon (2015) investigated the torques experienced by a 
circumstellar disc in an inclined circular binary, which is a 
special case of the general binary eccentricities studied in 
this paper. Nixon & Lubow (2015) demonstrated that an 
eccentric binary can exert non-zero Lindblad torques on a 
retrograde circumbinary disc, which are significantly weaker 
than those applied to a prograde disc, but can nonetheless 
be responsible for clearing a cavity around the binary. This 
is a special case of the arbitrary inclination framework used 
in this paper, although they did not make the assumption of 
a Keplerian disc. We do not specifically address retrograde 
discs in this paper. 

The most important technical aspect of our work is the 
computation of the Fourier components of the binary po¬ 
tential for arbitrary inclinations. We represent the poten¬ 
tial semi-analytically as a power series in eccentricity (as 
in the appendix of AL94), evaluated to order This is 
accomplished through the use of Wigner matrices to ac¬ 
count for arbitrary inclination between the disc and binary 
planes, and through an expansion of the binary orbit us¬ 
ing the Hansen coefficients of celestial mechanics. We note 
that Nixon & Lubow (2015) utilized an exact, non-series 
approach for computing the gravitational potential for the 
special case of retrograde circumbinary discs. This method 
amounts to numerically evaluating the Hansen integrals. 

This paper is organized as follows. In Section [2] we de¬ 


compose the gravitational potential of a misaligned binary 
into Fourier components with various forcing frequencies. In 
Section [ 3 ] we review how these potential components exert 
Lindblad torques on a disc. In Section |4] we apply this the¬ 
ory to the truncation of circumstellar discs, determining the 
size of the disc outer edge. We then determine the inner 
radii of circumbinary discs in Section [5] We address the va¬ 
lidity of our Keplerian disc approximation in Section [6] We 
summarize and discuss our results in Section [ 7 ] 


2 POTENTIAL COMPONENTS 

We consider a binary consisting of masses Mi and M 2 with 
total mass Mtot = Mi + M 2 , mass ratio q = M 2 /Mtot, semi¬ 
major axis a and eccentricity e. The orbital frequency is 
Hb = [GMtot/cd)^^^■ The orbit is inclined relative to a 
reference plane, taken to be the plane of a thin circumstellar 
or circumbinary disc, by an angle i. We will express the 
disturbing potential (per unit mass) acting on the disc in 
the form 

4? = ^ (r) cos [m<j) - (/r -f n) Hat] 

m,u.,n , , 

( 1 ) 

= ^ 4?^,^ (r) cos (mcj) — Nflst ), 

m,N 

where (r, (p) specifies the radial and azimuthal position of the 
disc particle, m is the azimuthal number in the disc plane, 
/i is the azimuthal number in the binary plane, A is a time 
harmonic number, and n is related to the eccentricity depen¬ 
dence of each potential component (see later in this section). 
In our formulation, m > 0, but and n can be positive, 
negative, or zero. Since different (/r, n)-components with the 
same (j) and t dependence are not physically distinct, we have 
defined 

4?m,iV (^) ~ ^ ^ ^it+n,N^7n,ii,n ix ) ; (2) 

which measures the total strength of the potential compo¬ 
nent having the azimuthal number m and rotating with the 
pattern frequency 

N 

cup = —Hb. (3) 

m 

Each component, denoted by (m,N), excites density waves 
at resonant locations in the disc, which give rise to the 
torques, Tm,N- We use different approaches to compute the 
potential ^rn,N, depending on whether we are considering a 
circumstellar or circumbinary disc. 


2.1 Circumstellar Disc 


We work in the reference frame centered on the primary 
star Ml, so that the disc rotation rate around Mi is Q (r) — 
(GMi/RY^^■ The gravitational potential felt by a disc par¬ 
ticle due to M 2 is 


$ = - 


GM2 

Ir-ral 


r • To 

+ GM2^^, 

^12 


( 4 ) 


where r is the particle’s position vector, and r 2 = ri 2 r 2 is 
the position vector of M 2 relative to Mi. The second term 
in the potential is the indirect part arising from the motion 
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of Ml around the center of mass of the system (e.g., Mur¬ 
ray & Dermott 1999). The direct term can be expanded in 
Legendre polynomials, leading to 

OO I 

-l> = -GM2^^P,(f.f2). (5) 

1=2 ’’12 

The Z = 0 term is an irrelevant constant which has been 
dropped. The angular dependence can be separated by an 
expansion in spherical harmonics, 

i 

■ ^2) = E {02, h) Yl,m {0, D ■ (6) 


We define two coordinate systems. The unprimed coordinate 
system (6, (j>) has its z-axis aligned with the disc angular 
momentum, z = Ld. The primed coordinate system {6', (j>') 
has its 2 :-axis aligned with the binary angular momentum, 
z' = Lb. The two coordinate systems share a j/-axis, and z • 
z' = cos (i). The spherical harmonics in these two coordinate 
systems are related through the Wigner d matrices. 


i 

Yi,^ ( 62 , ^ 2 ) = E 4,- (*) Yi,^ (Of ^ 2 ) ■ (7) 

fj,= — I 


Combining equations (O, di) and 0, evaluating Yi^rn in the 
disc plane {6 = 7r/2) and in the binary plane (O' = 
7r/2), and taking only the real part, the potential ((S]) can be 
written as 


E E E (i) 


where 


1 = 2 771=1 fj. = — l 

X (0 (“) ^ ^ cos {m(j} - fi(j> 2 ) , 

(Z — m)!l 


Wl,m = 


pr ( 0 ) 


_(Z -I- m)!_ 

= [(Z - m)! (Z + m)!] 2 

-1 


I fl-m 


I + m 


( 8 ) 


( 9 ) 


with the factor (—taken to be zero if Z -|- m is odd, 
and similarly for The product is zero unless 

Z, m and are all even or all odd. Notice that we have 
taken twice the sum over only positive values of m, since 
the (—m, —fi) and (m, fi) terms are identical (those with 
positive or negative /r are still distinct). In this form, the 
expression for the potential is exact, given an explicit forms 
of ri 2 /a and (fif the radial coordinate and true anomaly of 
M 2 in the binary plane. As a final step, we use their elliptic 
expansions to write (see Appendix 1X1) 



{m4> - = 


00 

E COS [m(j) - (n + n) ^bZ] ■ 

77 = —00 


( 10 ) 


Each coefficient is a series in powers of e, with the 

leading term proportional to The main approximation 
we make is to truncate these coefficients at a finite order 


in eccentricity (e*^®). The final expression for the potential 
strengths is 




2 GM 2 

a 




where /min = max (m, |/i|, 2). 



(11) 


2.2 Circumbinary Disc 


For circumbinary discs, we work in the reference frame cen¬ 
tered on the center of mass of the binary. In this frame, 
the orbital frequency of a particle in the disc is fl (r) — 
(GMtot/Dy^^. The disturbing potential can then be ex¬ 
pressed as 


<^ = -E 


GMi 


cfro"" 


Pi (f ■ f2) , 


( 12 ) 


where ri 2 is the separation between Mi and M 2 , and 

Ml = [g (1 - qf + (-1)' (1 - q) f] Mtot (13) 

(e.g.. Ford, Kozinsky & Rasio 2000; Harrington 1968). In 
general, odd-Z components are weaker than even-Z compo¬ 
nents for similar mass binaries, and for equal mass binaries 
{q = 1/2), Ml is identically zero for odd Z. Again writing 
the Legendre polynomials in terms of spherical harmonics 
in the binary and disc frames using the Wigner functions 
[equations ([51) and (0], we have 


$ = _2^ ^ ^ ^ (i) 


a 

1 = 2 771=1 ^= — 1 

' r \ “C^+i) / ri 2 \ ^ 


(14) 


^ (a) (^) ■ 


We expand the time-dependent orbital coordinates of the 
binary in a manner analogous to the circumstellar disc case. 


[ —cos [m(j) - = 


\ a / 


E cos [m4> - (n + n) HbZ] , 


(15) 


so the potential component is given by 

GM,,, ^CB ,i 




See Appendix[^for the values of the coefficients. 


(16) 


3 LINDBLAD TORQUES 

Each potential component $m,]v, rotating with the pattern 
frequency tup = NQ.^/m, excites density waves at the Lind- 
blad resonances (LRs), where 

tup - fZ (r) = (17) 

m 

where the upper (lower) sign corresponds to the outer (in¬ 
ner) LR. From here on we assume that the epicyclic fre¬ 
quency K (r) is equal to fZ (r) and both are proportional to 
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r i.e., the disc is exactly Keplerian. The LRs correspond 
to locations where 


and are located at 

t'lr J i 1)(1 — 5 )^^^ circumstellar disc 

ffl 1 [(m ± 1 ) circumbinary disc 

(19) 

The torque on the disc at a LR is (Goldreich & Tremaine 
1979) 


Hb 


N 


mil’ 


(18) 






r^LR, 

^ m,N 


dP 

dlnr 


1 ^. 


( 20 ) 


where E is the disc surface density, D = — rrd (D, — cjp)^, 

and 


^m,N 


d$7Ti,iv ^ 211 

d In r fl — cjp 


( 21 ) 


We also define which is the same as the above ex¬ 
pression but with replacing ’hm.iv- While ®m,jv is 

the quantity which determines the torque, it is useful to dis¬ 
cuss how various components contribute to it, since 

these are physically associated with different azimuthal forc¬ 
ing components of the binary orbit. Note that since we as¬ 
sume Keplerian discs, 211/ (H — cjp) = ^2m at r = tlr, 
and 


dP 

dlnr 


= T 


3N^ 

mil 


nl. 


( 22 ) 


At inner Lindblad resonances (ILRs), < 0, i.e., disc 

particles lose angular momentum, and at outer Lindblad 
resonances (OLRs), > 0, so disc particles gain angular 
momentum. Torques are also applied to the disc at coro¬ 
tation resonances, where wp — H (r) = 0 , however these are 
not important in disc truncation. Therefore we subsequently 
drop the superscript “LR” on Tm,N, as we only consider 
Lindblad torques. 

The viscous torque on the disc, assuming the a-ansatz 
for the kinematic viscosity coefficient, v = aCg/ll, is given 
by (e.g., Pringle 1981) 


(23) 


where h = H/r is the disc aspect ratio. As in AL94, we 
assume a gap is opened at the {m,N) LR if \Tm,N\ |2T|. 
Throughout this paper, unless otherwise noted, we adopt a 
disc model with h = 0.05 and a = 0.01. 


4 RESULTS: CIRCUMSTELLAR DISC 

4.1 Resonances Relevant to Outer Disc 
Truncation 

We are interested in the largest outer radius rout that a cir¬ 
cumstellar disc can have, given the orbital parameters (a, 
e, and i) of the binary. This amounts to determining the 
smallest radius at which a particular resonant torque can 
open a gap. Therefore we should consider the inner Lind¬ 
blad resonances of the potential components with the largest 
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Figure 1. The potential component ’l>m,N = 'I' 2,2 (solid line) 
and the dominant sub-components which contribute to it 

(dashed lines) for a circumstellar disc in an equal mass binary 
with e = 0.5. 


pattern frequencies that satisfy iTu.ivI ^ \Tv\, with |T„| eval¬ 
uated at the resonance location. The pattern frequency as¬ 
sociated with "hm.jv is cup = NQ.p/m, so we should examine 
components with small m and positive N. The smallest m- 
components whose ILRs are in the disc have m = 2, since all 
m = 1 ILRs are formally located at the origin [see equation 
m]- Therefore we focus on the (m, N) — (2, N) ILRs with 
N ^ 2, which are located at the fl/flB = N-.l commensura- 
bilities. 


4.2 Effects of Disc Inclination 

An example of the different contributions to a potential 
component is depicted in Figure [T] Here = 4^2,2 is 

shown (solid line) for an equal mass binary with e = 0.5 as 
a function of disc-binary inclination i. The potential is the 
sum of all 4'2,;i,n’s with /r -|- n = 2. The dashed lines de¬ 
pict each of the individual components which significantly 
contribute to 4^2,2 ■ Other components, not shown in Figure 
[T] such as V[' 2 ,- 4 , 6 , contribute negligibly (for e = 0.5), but 
can be important at larger e. We include components up to 
jnj = 8 in all subsequent calculations. At small inclinations, 
4^2,2 ~ 4/2, 2 , 0 , an equality which is exactly true at i = 0°, 
since (0) = 5^^^- In other words, at small inclinations, 
the = 2 component of the binary potential is primarily 
responsible for exciting the m = 2 disturbances in the disc. 
Also note that this is the only non-zero component when 
e = 0, since the others have jnj > 0, and oc At 

intermediate inclinations, the azimuthal components of the 
binary orbit other than /r = 2 become important, especially 
the /i = 0 term (4/2, 0 , 2 ), which contributes strongly between 
about 30° and 150°. At inclinations close to 180°, the com¬ 
ponent 412 ,- 2,4 becomes dominant. This component, despite 
being somewhat suppressed by its eccentricity dependence, 
becomes significant due to the strong coupling of the /r = —2 
and m = 2 terms near counter-alignment. Thus, an eccentric 
counter-aligned binary can exert non-zero Lindblad torques 
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Figure 2. Ratio of resonant torque to viscous torque as a function of inclination for the (m, N) = (2, N) inner Lindblad resonances 
for a circumstellar disc in an equal-mass binary. Each panel corresponds to a different binary eccentricity. The horizontal dashed line 
corresponds to |T^ jvl = Tv\, above which a given resonance can clear a gap. The dashed red lines (in the middle and right panels) are 
the (2, 2) torques computed using only the 'I'2,2,o potential component, rather than a sum of all appropriate 'I'm.fi.n’s. 


on a circumstellar disc. This is also true for circumbinary 
discs (see Section [S]). 

The inclination dependence of each component is deter¬ 
mined by a weighted sum of Wigner functions, but to leading 
order is determined by the first term in this sum. For exam¬ 
ple, ® 2 , 2 ,o oc (*) «: cos"‘ (i/2) and 'If 2 ,o ,2 oc do ,2 (*) oc 
cos^ (i/2) sin^ (i/2), both dominated by the I = 2 terms. 
However, the leading order term is not necessarily a good 
approximation to the actual i dependence. For example, the 
leading term of v!' 2 , 2 ,o is reduced (compared to i = 0°) by 1/4 
at 90°, while including all terms shows that this reduction 
is actually about 1/10, a much steeper decrease with incli¬ 
nation than given by d 2,2 (*)• The shape of the total ' 1 ' 2,2 is 
therefore a sum of many different Wigner functions, leading 
to its characteristic shape. There are two unique characteris¬ 
tics of this curve. First, it changes sign twice, at i = 29° and 
i = 172° (the particular values of these angles depend on 
the binary eccentricity), so that no torque is exerted on the 
disc at these inclinations. The sign of '^ 2,2 is irrelevant since 
the torque is proportional to 14/2,21^. Second, properly sum¬ 
ming the relevant 'Ifm,,n,n’s can counterintuitively lead to a 
H'm.AT whose absolute value (and hence resultant torque) is 
smaller than (the only non-zero component for 

an aligned disc), as it does for a large range of i In Figure [T] 

Figure [2] shows how the resonant torques responsible 
for disc truncation (normalized by viscous torque) vary with 
inclination, for an equal mass binary with several different 
eccentricities. Notably, even for a circular binary (the left 
panel of Figure [T]), there are non-zero torques with N > 2 
when the disc is not aligned. For example, |T 2 , 4 | depends on 
4^2,4, which is equal to '^ 2 , 4,0 for a circular binary. To leading 
order this depends on i as ^ 4,2 (*) oc cos® (i/ 2 ) sin^ (i/ 2 ), 
which is largest at i = 60°, and zero at 0° and 180°. In our 
canonical disc model (with h = 0.05 and a = 0.01), |r 2 , 4 | 
is never strong enough to clear a gap when e = 0, but if a 
or h of the disc were slightly smaller, it would be possible 
for a certain range of inclinations. This is impossible when 
the disc and binary are aligned. The largest torque for a 


circular binary is |72,2|, which depends only on 4^2, 2 , 0 , a 
monotonically decreasing function of i (see Figure [T|). This 
torque is very large at small inclination (exceeding \Tv\ by a 
factor of over 10^), and gets weaker as inclination increases, 
becoming about 18 times smaller (compared to aligned) at 
90°, and falling to zero at 180°. 

For e > 0 (the middle and right panels of Figure [2]), 
higher N torques generally become stronger (compared to 
e = 0). At small inclinations, 172, 2 I > |72,3| > |72,4| and so 
on. This is because at small i these torques are primarily a 
result of the '^ 2 , 2 ,n -2 potential components, which have the 
approximately the same inclination dependence [d |,2 (9 ^o 
leading order], but different eccentricity dependence (e^“^), 
the latter of which determines their relative strengths. At 
larger inclinations, the contributions of the yf 2 potential 
components cause the torques to have more oscillatory in¬ 
clination dependences (as demonstrated in Figure [T]). The 
dashed lines in the middle and right panels of Figure[2]show 
|72,2| computed using only the ’F 2 , 2 ,o component, demon¬ 
strating the relative importance of these couplings. As a 
result of this behavior, at large i (e.g., above about 120° for 
e = 0.3), the torque |72,Ar| is no longer monotonic in N. 


4.3 Location of Outer Disc Edge 

The following procedure is used to compute rout, the loca¬ 
tion of the outer edge of a circumstellar disc, for a given 
disc-binary inclination i. For each resonance (ILR), labeled 
by (m, N), or by f2(riLR)/flB = N/{m — 1), we first compute 
the torque jTm.jvl as a function of e (see Figure[3]). Then we 
find the range of e for which \Tm,N\ > \T\ (so that a gap 
can be opened), for each {m,N). Then at every value of e, 
the outer radius of the disc is identified as the location of 
the gap-opening resonance located at the smallest radius. 

The torques \Tm,N\ (normalized by viscous torque at 
the ILR, \Tv\) for {m, N) = (2,4) and (2,5) as a function 
of e are shown in Figure |3] For small e and i, |Tm,jv| is a 
monotonically increasing function of e, so that finding the 
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Figure 3. Ratio of resonant torque to viscous torque as a function of e for (m, N) = (2, 4) and (2, 5) ILRs in a circumstellar disc, for 
several inclinations. The horizontal dashed line corresponds to |T^ ;v| = 




o 



3:1 

7:2 

4:1 

5:1 

6:1 

7:1 

8:1 


Figure 4. Size of a circumstellar disc (location of outer edge rout) in an equal mass binary as a function of eccentricity for several 
disc-binary inclinations. The disc parameters are h = H/r = 0.05 and a = 0.01. The ratios r2(r/i,fl)/f2B for the resonances considered 
are labeled on the right. The filled points indicate the eccentricities and resonance locations at which the outer disc radius changes, and 
the lines connecting them indicate the behavior of rout between these points. The points represented by crosses indicate the minimum 
eccentricities required to open a gap at the Q/Q^ = 7:2 commensurability (for inclinations matched to the lines by color). This is an 
example of a resonance which is not important for disc truncation, since other resonances can open gaps at smaller radii for the same 
value of e. From left to right (in each panel), the four dashed lines correspond to: the size of a particle disc determined by Pichardo, 
Sparke & Aguilar (2005) (labeled “PSA05”), the long-term stability limit for S-type planets in binaries (for i = 0°) from Holman & 
Wiegert (1999) (labeled “HW99”), the average Roche lobe radius evaluated at the pericenter separation of the binary [labeled “r£i(e)”, 
see equation ||24| |1. and the binary pericenter separation, r 2 ,inin = (1 ~ which sets a strict upper limit for the disc size. 


range for which a gap can be opened amounts to finding 
the minimum e for which \Tm,N\ > \Tv\. However, for large 
inclinations, [Tm.jv| has local extrema which result in the 
existence of multiple values of e for which |Tm,iv| = \Tv\- For 
example, at t = 45° and 90°, as e is increased, \Tm,N/Tu\ 
exceeds unity above a critical e, then reaches a maximum, 
decreases, and again drops below unity at a second critical 
e. Another minimum is then reached, and then it exceeds 
unity again at yet another critical e. Despite the oscillatory 
features in Figure [3] once the (m, N) resonance is cleared, 
the {m,N -\- 1) resonance is always cleared before [Tu^^jvl 


becomes non-gap-opening. In other words, as e is increased, 
the innermost gap clearing resonance shifts inward, so that 
the disc size is a monotonically decreasing function of e for 
all inclinations we have considered. 

Figure |4] shows an example of the disc outer radius as 
a function of eccentricity for an equal mass binary, for sev¬ 
eral inclinations. In this figure, the filled points represent 
the eccentricities and resonant locations at which the outer 
radius abruptly changes due to a new innermost resonance 
“turning on” (becoming able to clear a gap). These result in 
the outer disc radius following the stairstep-shaped curves 
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Figure 5. Same as Figure [4] except for for discs with viscosity 
parameter a = 10“® and 10“^. The straight lines connecting the 
points are for graphical convenience only, and are a proxy for the 
starstep-shaped curves depicted in Figure [4] 



which connect the filled points. Several limiting radii are 
also shown in Figure (4] The dashed line labeled ri:,i(e) is 
the average radius of the Roche lobe around Mi (Eggleton 
1983) evaluated at the pericenter of the binary orbit: 


rLi{e) 


0.49gi^^(l-e)a 
0.6gi/®+ln + 


(24) 


where qi = Mi/M 2 = (1 — q)/q (with q = M 2 /Mtot). Equa¬ 
tion (|24ll is approximate since the original Roche lobe is 
calculated for circular, synchronized binaries. We see from 
Figure |4] that for i > 90°, our predicted disc radius (based 
on gap opening criterion) is larger than ri:,i(e), suggesting 
that our result for rout may be larger than the actual disc 
size. Note, however, that since equation dSH) approximates 
the Roche lobe as a sphere, it is possible for the disc to ex¬ 
tend beyond r_Li(e) for some inclinations. The dashed line 
labeled r 2 ,min corresponds to the pericenter separation of 
the binary [r 2 ,min = (1 — e)a]. This is a strict upper limit on 
the disc size when it is either aligned or counter-aligned (the 
disc-projected closest approach distance is slightly different 
when there is misalignment), since if one of the stars were 
allowed to plunge through the disc it would rapidly clear 
material from its orbit. 

In Figure 3] we have also shown the long-term stability 
limit for S-type planets in binaries (for i = 0°), based on the 
fitting formula of Holman & Wiegert (1999). This limit is 
smaller than the truncation radius of the disc, but close to 
the “orbit-crossing” limit determined by Pichardo, Sparke 
& Aguilar (2005). 

Our canonical disc model has h = 0.05 and a = 0.01. 
The effect of varying the disc properties is shown in Figure 
[5l which depicts the disc outer radius as a function of ec¬ 
centricity for an equal mass binary (as in Figure |4]) when 
a = 10“® and a = lO”'^. These correspond to a reduction of 
all viscous torques Ti, by a factor of 10 and 100 compared 
to our canonical model. Since is proportional to aD, the 


curves labeled a = 10“® represent any disc with a combi¬ 
nation of a and h such that (a/10~^) (/i/0.05)^ = 1, and 
similarly for curves labeled a = Although the viscous 

torques rescale trivially in this way, the complicated depen¬ 
dence of the resonant torques on eccentricity and inclination 
does not allow a simple rescaling of the disc size versus ec¬ 
centricity relations. This explains the qualitative differences 
between the curves in Figures |4] and [5] Most notably, for 
the values of a in Figure O a gap can be opened at the 
4:1 commensurability at zero eccentricity for i = 45° and 
90°, contrary to our canonical disc model (see Figure |4]). 
However, broadly speaking, the difference in disc size due to 
an order of magnitude change in aD is comparable to the 
difference in size due to a change of inclination of about 45°. 

The effect of binary mass ratio on circumstellar disc size 
is explored in Figure |6l Recall our definition of the mass 
ratio q = M 2 /Mtot, and that M 2 is always considered the 
perturber, regardless of whether it is more or less massive 
than Ml. Therefore the disc sizes for g = 0.1 and g = 0.9 
can be thought of as the sizes of the circumprimary and 
circumsecondary discs in a system in which the secondary is 
1/9 the mass of the primary, and similarly for g = 0.3 and 
g = 0.7 (in which case the secondary is 3/7 the mass of the 
primary). The behavior is qualitatively similar to the equal 
mass case (see Figure |4]), but with more massive perturbers 
leading to smaller discs and vice versa. 

We have restricted our calculation of rout to binaries 
with e < 0.8. At high eccentricity, our ability to compute the 
disc outer radius is restricted by our finite expansion of the 
disturbing potential: the highest order resonance we consider 
is H/Hb = 8:1. For larger values of e, higher order resonances 
can clear gaps, further reducing the size of the disc, but 
computing the appropriate potential components for these 
resonances is impractical in our semi-analytic approach. 


5 RESULTS: CIRCUMBINARY DISC 

5.1 Resonances Relevant to Inner Disc Truncation 

Determining the size rin of the inner cavity of a circumbi¬ 
nary disc amounts to finding the largest radius at which 
a gap can be cleared. Therefore we consider the OLRs of 
the potential components with the smallest possible pattern 
frequencies. Since H(roLR)/HB = N/{m + 1) [see equation 
m], we shall focus on A = 1 and m ^ 1, corresponding 
to the H/Hb = l:(m-|- 1) commensurabilities. However, for 
equal-mass or nearly-equal mass binaries, the strength of the 
odd-m potential components is zero or small relative to the 
even-m components, resulting in a relatively weak torque 
|7m_jv|. As we will also show, [Ti^ij can be weak even for bi¬ 
naries which are not close to equal mass. Therefore, it is also 
necessary to consider the next strong resonance interior to 
the 1:2 commensurability, namely the 2:3 commensurability 
(m = 2, N = 2), which can be cleared by the T 2,2 torque. 


5.2 Effects of Disc Inclination 

Figure 0 shows the resonant torques normalized by vis¬ 
cous torque for a circumbinary disc around a binary with 
mass ratio g = 0.3 (cf. Figure [2] for a circumstellar disc). 
The {m,N) = (1,1) torque is weak due to the absence 
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Figure 6. Same as Figure[4l except for several different mass ratios. As in Figure[5l the straight lines connecting the points are only for 
graphical convenience, and should be interpreted as stairstep-shaped depicted in Figure [4| 
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Figure 7. Ratio of resonant torque to viscous torque ratio as a function of inclination for the (m, N) = (m, 1) and (2, 2) outer Lindblad 
resonances in a circumbinary disc with mass ratio q = 0.3. Each panel corresponds to a different binary eccentricity. As in Figure!^ the 
horizontal dashed line corresponds to = l^i^l (above which a resonance can clear a gap), and the dashed lines in the middle and 

right panels are the (2, 2) torques computed using only the ^^2,2,0 potential component. 
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Figure 8. Ratio of resonant torque to viscous torque as a function of e for (m, N) = (2,1) and (3,1) outer Lindblad resonances in a 
circumbinary disc, for a mass ratio q = 0.3. The dashed horizontal line indicates \Tm,N\ = |Th|. 


of the dipole term in the binary potential [see equation 
mi)], so that to leading order, is proportional to 

dl l (*) Meanwhile '52,2 is approximately propor¬ 
tional to d% 2 (*) (and has its OLR at a smaller ra¬ 

dius), and so the (m,N) = (2,2) torque is generally much 
stronger than the (1,1) torque and can potentially be the 
most relevant to clearing the inner cavity at low e and i. The 
(2, 2) component is also relevant for an equal-mass binary, 
for which all torques with odd m are zero. 

For the same reasons as for a circumstellar disc, the 
torque can be an oscillatory function of i, but this effect is 
even more pronounced for a circumbinary disc. For example, 
for a circular binary, Ts,! can clear a gap at the 1:4 commen- 
surability for i between 37° to 106°, which is impossible for 
an aligned disc. For e > 0, there are several ranges of i for 
which the ordering of torques is very different than for an 
aligned disc, for which the torque strengths depend mono- 
tonically on the resonance location (innermost resonances 
are the strongest). 


5.3 Inner Cavity Size: Inner Disc Radius 

The procedure for computing the inner disc radius rin for a 
circumbinary disc at a given inclination i is as follows. First, 
for each resonance {m,N) = (m, 1) and (2,2), we compute 
|7m.jv| as a function of e (see Figure [S]), and determine the 
range of e values where gap opening is possible (|r,n,jv| ^ 
\Ttj\). The inner disc radius is determined by which gap¬ 
opening resonance is located at the largest radius for every 
value of e. 

As in the circumstellar disc case, the complicated de¬ 
pendence of \Tm,N\ on e can result in several ranges of e for 
which a resonance can open a gap. But unlike the circumstel¬ 
lar disc case, this effect is signihcant enough to qualitatively 
alter the dependence of the disc size on e. As an example, 
\Tm,N/Tu \ is shown in Figure[8]for (m, N) = (2,1) and (3,1), 
corresponding to the = 1:3 and 1:4 commensurabil- 

ities, for a binary with mass ratio q = 0.3, for several disc 
inclinations. While iTa.i] is a monotonically increasing func¬ 
tion of e, such that for each i there is a critical e above which 


it can open a gap, the behavior of |r 3 ,i| is very different. 
For i = 45° and 90°, it can open a gap at zero eccentricity, 
and it decreases with e so that above a certain critical e a 
gap can no longer be opened (for i = 45° a gap can again 
be opened at yet another larger e). In general, a given reso¬ 
nance, (m, N), is not always gap-opening at the lowest value 
of e for which the next resonance, (m + 1, N), first opens a 
gap. This is in contrast to the case of a circumstellar disc, 
for which gaps are always cleared at sequentially higher or¬ 
der resonances as e is increased. At large inclinations, this 
vastly different behavior of the jTm.jvl’s can result in cav¬ 
ity sizes which both increase or decrease with increasing e, 
sometimes exhibiting both behaviors for a single inclination. 

The inner disc radii for three different mass ratios and 
various inclinations are shown in Figure |5] We hrst focus on 
the middle row, which corresponds to the mass ratio q = 0.3 
(as in Figures [7] and [SJ. For f = 0° and 22.5°, the inner 
edge is located at the 1:2 commensurability for e = 0 and 
increases with e as gaps are opened at higher-order reso¬ 
nances. The maximum cavity size is located at the 1:6 com¬ 
mensurability i = 0°, and at the 1:5 commensurability for 
i = 22.5°. For i = 45° and 90°, the 1:4 commensurability can 
be cleared at zero eccentricity (see the left panel of Figure 
[71, but the cavity size then decreases with e as the 1:4 torque 
becomes too small and the 1:3 torque becomes responsible 
for gap opening. For i = 45°, the size then goes back up to 
the 1:4 commensurability at e = 0.47 and remains at that 
size for larger e. For i = 90° the cavity shrinks as e increases, 
then remains at the 1:2 commensurability. Clearly, the be¬ 
havior of the disc inner cavity size is more complicated than 
in the circumstellar disc case, and even its qualitative be¬ 
havior (increasing or decreasing with e) cannot be trivially 
ascertained, but requires a full consideration of the details 
of the resonant torques. 

The top and bottom rows of Figure O give the inner 
disc radii for g = 0.1 and g = 0.5. Note that a mass ratio g 
greater than 1/2 is degenerate with a mass ratio of 1 —g for a 
circumbinary disc. Much of the qualitative behavior seen in 
the g = 0.3 case can also be seen for g = 0.1, for example the 
cavity size is an increasing function of e for i ^ 45°, while 
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Figure 9. Size of inner cavity (location of inner edge rjn) of a disc around a binary with several mass ratios (corresponding to different 
rows), each for several binary-disc inclinations. The disc parameters are h = 0.05 and a = 0.01. The ratios fl(roLR)/^B foi" 
resonances considered are labeled on the right. The filled points indicate the eccentricities and resonance locations at which the inner 
disc radius changes, and the lines connecting them indicate the behavior between these points. In the middle row, the points represented 
by crosses indicate the minimum eccentricities for which a gap can be cleared at the 2:7 commensurability, for i = 0° and 22.5° (for 
higher inclinations, a gap cannot be cleared for any value of e). This is an example of a resonance which is not important in determining 
the size of the inner disc, since gaps can be opened at larger radii for the same values of e. The dashed lines in the left panels labeled 
= Q” indicate the radius below which the disc does not satisfy the Rayleigh stability criterion (for i = 0°). The dashed lines labeled 
‘V 2 ,max” correspond to the radial location of the secondary component of the binary at apocenter [r 2 ,max = (1 — g)(l + e)a], which is 
a strict lower limit on rjn (for q = 0.5, r 2 ,max is below the displayed range of the y-axis). The dashed lines labeled “HW99” indicate 
the stability limit for P-type planets in binaries (for i = 0°) from Holman Sz Wiegert (1999). Note that for q = 0.3 and q = 0.5, when 
2 = 0° and i = 22.5°, the cavity transitions from its smallest size (the 1:2 commensurability for q = 0.3 and the 2:3 commensurability 
for q = 0.5) to a larger size (1:3 in both cases) at a very small but non-zero eccentricity (e = 0.01 — 0.02). 


© 0000 RAS, MNRAS 000, 000-000 











































































Tidal Truncation of Inclined Circumstellar and Circumbinary Discs in Young Stellar Binaries 11 


for i = 90°, it can either be decreasing or increasing over 
different ranges of e. A new behavior is seen for i = 135° and 
q = 0.1: the disc size is indepedendent of e. In this case, for 
all values of e, the torque iTi^l is strong enough to clear the 
1:2 commensurability, but no other torque iTm.jvl can clear a 
gap at its resonant location in the disc regardless of the value 
of e. For an equal mass binary [q = 0.5), there are fewer 
possible truncation sites because odd m torques are zero 
(since Mi = 0 for odd 1). The only allowed cavity sizes are at 
the 2:3, 1:3 and 1:5 commensurabilities. The dependence of 
the cavity size on e is simple in this case: for all inclinations, 
Tin increases with e, with the maximum cavity size located 
at the 1:5 commensurability for i = 0° and i = 22.5° and at 
1:3 commensurability for other inclinations (45°, 90°, 135°). 

Figure [9] also shows, as a function of eccentricity, the 
radius at which the squared radial epicyclic frequency [see 
equation (I26II ] is equal to zero. Below this radius, disc par¬ 
ticles are unstable to radial perturbations (according to the 
Rayleigh stability criterion), and so the disc can not ex¬ 
tend inwards beyond this radius. This radius is shown in 
the left panels, and is evaluated for i = 0° (or equiva¬ 
lently, i = 180°). It is not shown for larger inclinations, 
for which the instability region shifts inward and lies inte¬ 
rior to r 2 ,max = (1 — (?)(! + e)o, the radial coordinate of 
the secondary star at the apocenter of the binary orbit (also 
shown in Figure [Oj, which is a strict lower limit for rin. For 
the mass ratios and inclinations we have considered, nei¬ 
ther of these non-resonant truncation mechanisms constrain 
the inner cavity size, except for at very large e for g = 0.1 
and i ^ 90° (see the two upper right panels), for which the 
Lindblad torques are weak. Note that for a retrograde disc 
{i = 180°), for which the torques are further weakened, the 
Rayleigh stability criterion may be more relevant for trun¬ 
cation than gap clearing by Lindblad torques for a wider 
range of binary parameters (see Nixon & Lubow 2015). 

The long-term stability limit for P-type planets in bina¬ 
ries (for i = 0°) (Holman & Wiegert 1999) is also shown in 
Figure [5] The inner truncation radius is inside the stability 
limit, which is of interest to the formation of observed cir¬ 
cumbinary planets, many of which reside very close to the 
stability limit (see Welsh et al. 2014 and references therein). 
Thus, planets that have formed near or migrated to the in¬ 
ner edge of the disc (e.g., Kley & Haghighipour 2014) may 
experience dynamical instability as the disc disappears. 


6 EFFECTS OF NON-KEPLERIAN ROTATION 

Up to this point, we have considered discs with Keplerian 
rotation profiles (fl oc for which the rotation fre¬ 

quency n(r) is equal to the radial epicyclic frequency K.{r). 
Strictly, this is valid when the static, axisymmetric part of 
the gravitational potential consists only of a monopole term, 
proportional to 1/r. However, in general, the rotation profile 
is modified by the quadrupole and higher-order components 
of the binary potential^ The squared rotation frequency is 


^ Nixon & Lubow (2015) included this effect in their analysis of 
the Lindblad torques experienced by a retrograde circumbinary 
disc. 


given by 

^\r) = -^[$K{r) + ^oAr)], (25) 

r dr 

and the squared radial epicyclic frequency by 

Here <l>K(r) is equal to —GMi/r for a circumstellar disc 
and —GMtot/r for a circumbinary disc, and $ 0 , 0 (f) is the 
(m, N) = (0,0) component of the appropriate disturbing po¬ 
tential. Equations (I25II and (1261) assume that pressure gradi¬ 
ents, which are an additional source of non-Keplerian rota¬ 
tion, contribute negligibly compared to gravity (see Section 
0 . 

The non-Keplerian rotation profile has two effects on 
the Lindblad resonances. First, the resonance locations are 
shifted, so that they no longer correspond to exact integer 
commensurabilities of t}{r) and fie [as in equation (1191) ]. In¬ 
stead, they are determined by solving the general LR condi¬ 
tion [equation (1171) 1. and are functions of e and i. We adopt 
the notation for the Keplerian LR locations [as given 
by equation m] to distinguish them from exact locations, 
ruR. Second, the Lindblad torques [equation (1201) ] are mod¬ 
ified from their values in a Keplerian disc, since dD/dlnr 
and ^rn,N must be evaluated at the new resonance loca¬ 
tions using the modified expressions for fl and k. We also 
adopt the notation jy for f ho torque on a Keplerian disc 
to distinguish it from the true torque In this section, 

we consider how these effects modify the properties of the 
resonances relevant to truncation for both circumstellar and 
circumbinary discs. 

6.1 Circumstellar Disc 

Figure [in] (left panel) gives an example of how the disc rota¬ 
tion profiles and the locations of ILRs (for m = 2) are mod¬ 
ified by a non-Keplerian potential. The profiles are obtained 
by evaluating equations (ESJ - (uni) numerically. At high ec¬ 
centricities, the strong deviations from Keplerian rotation 
render some resonances non-existent, since there is no r for 
which 17 —ft/m is equal to cup. This occurs for e > 0.4 for the 
N = 3 ILR, and for e > 0.6 for the N = 4 ILR. Notably, the 
(m, N) = (2, 2) ILR (nominally the 2:1 orbital commensura¬ 
bility) does not exist for an aligned disc in a circular, equal 
mass binary. The largest fractional shift in rpR for a given 
resonance is about 12%, occurring for the largest value of e 
for which the resonance still exists. For eccentricities near 
this value, there can also be a second location satisfying the 
resonant criterion. However, the second location is unlikely 
to be relevant to disc truncation, as it is close to the point 
at which the disc becomes Rayleigh unstable, and is located 
at a larger r than the original resonance. 

The numerical results for the ratios of resonant torque 
Tm,N to viscous torque T^, are depicted in Figure [TT] along 
with their Keplerian counterparts. Note that the correction 
to Ti, at the modified resonance location has been taken into 
account in Figure[Tl] The deviation of the torque ratio from 
its Keplerian counterpart is largest for i = 0°, and negli¬ 
gible for other inclinations (comparable to the thickness of 
the plotted lines). The i = 0° curves terminate at a smaller 
value of e than for other inclinations, due to the fact that 
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Figure 10. The resonant frequencies (combinations of rotation frequency Q and radial epicyclic frequency k) which determine the 
location of several Lindblad resonances [see equation (Il7l l1: Q — k/ 2 (in units of the binary orbital frequency ^b), which is relevant to 
m = 2 ILRs for a circumstellar disc (shown for an equal mass binary with i = 0°; left), and Q + k/ 2, which is relevant to m = 2 OLRs 
(shown for a circumbinary disc with q = 0.3, for i = 0° and i = 90°; right). The line colors correspond to different binary eccentricities, 
as indicated by the color bar on the right. The curves terminate (at large r for a circumstellar disc, and at small r for a circumbinary 
disc) at the point at which k? becomes negative. The dashed curves indicate the corresponding frequencies (Q di k/ 2) in a Keplerian 
disc. The horizontal dashed lines represent pattern frequencies wp = NQq/2 for several values of N. The intersection of the Q di k,/2 
curves with the horizontal lines correspond to the locations of the (m, N) = (2, N) LRs. Note that for a circumstellar disc, Q and k. do 
not strongly deviate from their Keplerian values at large inclinations, therefore only the 2 = 0° case is shown. 



Figure 11. Ratio of resonant torque to viscous torque Ti^, both evaluated at the true resonance location rpR (solid lines), and 

the ratio of these torques in a Keplerian disc, evaluated at the Keplerian resonance location (dashed lines), for a circumstellar disc 
in an equal mass binary. The {m, N) = (2,3) and (2,4) ILRs (left and right panels), which are the lowest-order resonances responsible 
for truncation at non-zero eccentricity, are each shown for several inclinations. The horizontal dashed line corresponds to the threshold 
for disc truncation (|T^ ,N/n\ = i)- 


the resonance ceases to exist (see Figure fTOl) . For eccentrici¬ 
ties close to this point, the torque deviates strongly from its 
Keplerian value (by over an order of magnitue), due to the 
rotation profile becoming highly distorted. However, these 
large deviations occur at much larger values of e than the 
ones for which the resonance first truncates the disc. The 
eccentricities for which iTm.ivI = llhl change by less than 
0.01 for the resonances shown in Figure [TT] Higher order 
resonances, located at smaller radii, are affected even less 
strongly. Therefore, the assumption of a Keplerian disc is a 
reasonable approximation for determining outer disc trun¬ 
cation. 

To quadrupole {I = 2) and order, explicit expressions 


for the disc rotation profiles and Lindblad torques can be 
obtained. The axisymmetric, time independent part of the 
disturbing potential can be approximated by 

$o,o(r)«-^/cs(e,i)(0', (27) 

where the dependence on the inclination and eccentricity of 
the binary is given by 

/cs(e,i) = i [3cos^(i) - l] (^1 + • (28) 
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Thus, to this order, 

n2/ \ ^ r>2 I 24 >o,o(»') j 2 / \ 
U (r) « SiK H- 2 - ^ Y) 


^2 , 8'ho,o(r) 

“k+ - -2 -• 

(29) 


The resonance condition for m = 2 ILRs, which are relevant 
to outer disc truncation, is then 

JIk (1 + 2ecs) ~ A'TJb, (30) 


where we have defined the dimensionless parameter 

qfcs{e,i) 


ecs 


^’0,0 


4iV2 


The resonances are therefore located at 




^LR ( 1 + 


(31) 


(32) 


and hence shift outward (inward) compared to their loca¬ 
tions in a Keplerian disc when i is less (greater) than 54.7°. 
We can use equations (EH) and 1321) to evaluate the Lind- 
blad torque Tm,N [equation (1201) ] explicitly (to quadrupole 
and order). For m — 2, we find 

T2,n (r-LR) « (’'lr) ^1 + ^*=03^ . (33) 


6.2 Circumbinary Disc 


In Figure [10] (middle and right panels), we show the loca¬ 
tions of the (m, N) = (2,1) and (2, 2) OLRs for i = 0° and 
90°, including the effects of the non-Keplerian potential. The 
largest fractional shifts are about 4%, provided that the res¬ 
onance exists. The (m, N) = (2, 2) resonance does not exist 
for e > 0.5 when i = 0°, and for e > 0.7 when i = 90°. 

The ratio of resonant to viscous torques is shown in 
Figure [T^ The {m,N) = (1,1) torque differs most from the 
Keplerian case for i = 0°, especially near e « 0.6, above 
which the resonance no longer exists. For i = 90°, there is 
also some appreciable deviation from the Keplerian torque at 
large values of e. However, the values of e for which the res¬ 
onances become gap-opening change negligibly. Therefore, 
as for the case of a circumstellar disc, the approximation of 
a Keplerian disc does not affect truncation. 

To leading order in I and e, the non-Keplerian part of 
the potential can be approximated by 


■1*0,0 «-^^/cB(e, i) (-) , (34) 


with 


fcB{e,i) 



-1^2 • 2 
loe sin 



(35) 

where fj,B = MiM 2 /Mtot is the reduced mass of the binary. 
The resulting rotation frequency and epicyclic frequency are 




^2 3‘I>o,o(’') J 2 / N 
Hk - ’ and k (r) 


o2 , 3'l>o,o(r) 

H--• 


(36) 


The resonance condition for N = 1 OLRs, which are relevant 
to inner disc truncation, is 


Hk 



m — 1 
m -I- 1 


ecB 


Hb 

m + 1 ’ 


(37) 


where we have defined the dimensionless parameter 


ecB 




pK 

Tr 


So the resonance locations are 


g(l - <?)/cB(e,f) 

4(m + 1)^/® 


ruR 


K 

' >'lr 


, , m — 1 

1 + -— I ecB 

m -I- 1 


(38) 


(39) 


In this approximation, the {m,N) = (1,1) resonance (nom¬ 
inally the H/Hb = 1:2 commensurability) is not shifted rel¬ 
ative to its position in a Keplerian disc. The other = 1 
resonances shift to larger (smaller) radii compared to their 
locations in a Keplerian disc when /cB(e, i) is positive (neg¬ 
ative). Explicitly evaluating the Lindblad torque to the same 
order of approximation, we find that for N = 1 (and m > 1), 


Tr.,1 (ruR) « T^l (r^^a) 

m -I- 2 12m 

m -I- 1 3m -I- 1 


xil + 


2{m - 1) 


ecB 


(40) 


7 SUMMARY AND DISCUSSION 
7.1 Main Results 

We have developed a method for computing Lindblad 
torques due to a binary potential on circumstellar and cir¬ 
cumbinary discs which are misaligned with the binary or¬ 
bital plane. We used this theory to determine the outer radii 
of circumstellar discs and the inner radii of circumbinary 
discs, generalizing the work of Artymowicz & Lubow (1994; 
AL94) for aligned discs. The summary of our results is as 
follows. 

In the presence of misalignment (and non-zero eccen¬ 
tricity), each azimuthal component of the binary potential 
experienced by the disc is a result of many azimuthal compo¬ 
nents of the potential in the binary plane. This is in contrast 
to aligned disc, for which each azimuthal component of the 
potential can only produce perturbations in the disc which 
have the same azimuthal number. As a result, the inclination 
and eccentricity dependence of the gap-opening Lindblad 
torques can be somewhat complicated. Rather than sim¬ 
ply decreasing with inclination i (as for circular binaries), 
or increasing with eccentricity e (as for aligned discs), the 
Lindblad torque, |r,u,jv|, associated with the potential com¬ 
ponent which rotates with pattern frequency ujp = NQb Im 
(where Hb is the binary orbital frequency, and m, N are 
integers), generally has multiple extrema as functions of e 
and i. 

As in AL94, we adopt the resonance gap opening crite¬ 
rion \Tm,N\ > \Tv\, where is the viscous torque evaluated 
at the resonance location. For circumstellar discs, the most 
important resonances for disc truncation are located at the 
H/Hb = A:1 commensurabilities, with N ^ 2. Ignoring non¬ 
resonant truncation mechanisms, the outer edge of the disc 
is determined by the innermost of these resonances which 
can clear a gap. Despite the complicated dependence of in¬ 
ner Lindblad torques |Tm,jv| on inclination and eccentricity, 
for a given inclination, the resultant outer disc radius is a 
decreasing function of eccentricity (see Figures |4}]6]|. This is 
the same qualitative behavior as for an aligned disc (AL94). 
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Figure 12. Same as Figure [TT] but for the (m, N) = (1,1) and (2,1) OLRs in a circumbinary disc with mass ratio q = 0.3. 


Larger inclinations lead to larger discs: in an eqnal mass bi¬ 
nary, the disc is about 20% larger for i = 90° and 40% larger 
for i = 135°, compared to an aligned (i = 0°) disc. 

If the innermost gap-opening resonance lies outside of 
the tidal radius, then the disc is truncated non-resonantly. 
This tidal radius has been estimated to be 75 — 90% of 
the Roche lobe radius for aligned discs (Paczyfiski 1977; 
Papaloizou & Pringle 1977; Paczyhski & Rudak 1980; 
Pichardo, Sparke & Aguilar 2005). In the absence of an 
equivalent theory for misaligned discs, we estimate the tidal 
radius to be r_Li(e), the average Roche lobe radius at the bi¬ 
nary pericenter separation [see equation (I24ll 1. We find that 
a circumstellar disc in an equal mass binary fills its Roche 
lobe if inclined more than 45°, for our standard disc parame¬ 
ters {h = 0.05 and a = 0.01). If ah^ is smaller by a factor of 
10 ^ —10^ (a reasonable range for protoplanetary discs), reso¬ 
nant truncation can restrict the disc size to less than r_Li(e) 
for inclinations up to 90°. In non-equal mass binaries (we 
have considered M 2 /M 1 = 2.3 — 9.0), for our standard disc 
parameters, the circumprimary disc fills its Roche lobe if 
inclined by more than 45°, while the circumsecondary disc 
barely fills its Roche lobe at an inclination of 90°. Since 
the Roche lobe formula [equation (1241) ] is approximate, 3D 
numerical calculations are needed to determine the precise 
outer radii of such discs. 

We have also considered the resonance clearing of the 
inner cavity of a circumbinary disc. The most important 
resonances are the D/Db = 1:A^ {N ^ 2) and 2:3 commen- 
surabilities. As for a circumstellar disc, the outer Lindblad 
torques have a complicated dependence on inclination and 
eccentricity. However, unlike for a circumstellar disc, the de¬ 
pendence of the inner disc truncation radius on eccentricity 
can be very different at large inclinations (see Figure[9]). For 
an aligned disc, the inner cavity radius increases with eccen¬ 
tricity, and cannot extend past the 1:6 commensurability. 
This is also true at small inclinations (e.g., 22.5°), for which 
the cavity is slightly smaller than for an aligned disc. At 
larger inclinations, the cavity size can either be an increas¬ 
ing or decreasing function of eccentricity, possibly exhibiting 
both behaviors of different ranges of e. Nonetheless, the in¬ 
ner disc radius is generally smaller at large inclinations, for 
example, never extending past the 1:4 commensurability for 


i = 45° or i = 90°, or past the 1:3 commensurability for 
i = 135°. 

7.2 Approximations and Uncertainties 

The calculations presented in this paper adopt several as¬ 
sumptions. We have assumed the disc to be flat. This is 
reasonable, as under typical conditions, bending waves and 
viscous stress lead to small disc warps (see Foucart & Lai 
2014). Our results for the disc truncation radii (for both 
circumstellar and circumbinary discs) are based on the as¬ 
sumption of exactly Keplerian rotation profiles (D oc 
The validity of this approximation is addressed in Section 
[5] We find that, although non-Keplerian effects due to the 
binary potential can modify the resonance properties sig¬ 
nificantly at high eccentricities (e.g., some resonances can 
be rendered non-existent), gap-clearing by the relevant res¬ 
onances occurs at low eccentricities and is only slightly mod¬ 
ified. Thus, non-Keplerian rotation effects have a negligible 
effect on disc truncation. We note, however, that we did 
not consider the effect of strong pressure gradients near the 
edge of the disc, which may cause further modification of 
the rotation profile and LRs (Petrovich & Rafikov 2012). 

Our most important assumption (which is also made 
in AL94) is that angular momentum is deposited into the 
bulk disc material at, or very close to, the Lindblad reso¬ 
nances. This may affect our results in an appreciable way. 
First, angular momentum cannot be directly imparted to 
the disc material at Lindblad resonances, but must be car¬ 
ried as waves (Goldreich & Nicholson 1989), and received 
by the disc where the waves dissipate, either due to viscous 
damping (Takeuchi, Miyama & Lin 1996) or due to wave 
steepening and shock formation, although the latter should 
occur almost immediately after the waves are excited for the 
mass ratios we have considered (Goodman & Rafikov 2001). 
Second, the excited waves have broad angular momentum 
flux profiles, with widths comparable to their radial wave¬ 
lengths (e.g., Meyer-Vernet & Sicardy 1987), rather than 
being sharply peaked at the Lindblad resonances. For these 
reasons, the locations of disc truncation computed in this 
work may differ from those produced by self-consistent nu¬ 
merical treatments which account for these effects. For ex- 
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ample, the hydrodynamic simulations of discs aligned with 
equal mass circular binaries by MacFadyen Sz Milosavljevic 
(2008) show that while the resonance responsible for clear¬ 
ing the inner cavity (m = 2^ N = 2) is formally located at 
r = 1.31a for a Keplerian disc, the actual cavity extends to 
nearly 2a. This suggests that the details of wave excitation, 
propagation and breaking (or damping) are important un¬ 
certainties in the disc truncation process. With this is mind, 
our results for the truncation radii as functions of binary 
eccentricity and inclination should be interpreted as gen¬ 
eral trends rather than exact, sharply defined boundaries. 
We emphasize the need for detailed numerical work to fully 
interpret and assess the results of this paper. 
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APPENDIX A: EXPANSION OF POTENTIAL 
AS POWER SERIES IN ECCENTRICITY 

We expand (ri 2 /a)~^^^^^ cos (mt/) — ^,(^ 2 ) [foi" ^ circumstel¬ 
lar disc, see equation GQl)], or {r\ 2 /af cos {mrl) — u4>2) [for 
a circumbinary disc, see equation m] , as a power series 
in e. Here ri 2 is the binary separation and (f >2 specifies the 
azimuthal position of M 2 relative to Mi in the orbital plane 
of the binary, which, without loss of generality, can be cho¬ 
sen to be the true anomaly of the binary. We use the ellip¬ 
tic expansions (e.g., Murray & Dermott 1999; Brouwer & 
Clemence 1961) of rvija and ()) 2 , in terms of e and mean 
anomaly M = For the binary separation ri 2 we have 

1 °° 1 d 

— = 1 -I- -e^ - 2ey^ —— Js (se) cos (sM) 
a 2 ^s^de 

S = 1 

= 1 — e cos (M) + ^ [1 — cos (2M)] -I- [cos (M) — cos (3M)] 
2 8 
4 

+ — [cos (2M) — cos (4M)] -I- O (e®) . 

(Al) 

The following expression relating the mean anomaly M and 
eccentric anomaly E is useful: 


00 - 

E = M + 2J2-T (se) sin (sM) 


= M + e sin (M) -I- ^ sin (2M) -|- e® 


|sin(3M) 


+ e^ 


- sin (4M) - - sin (2M) 
3 6 


+ O (e®) . 


isin(M) 

(A2) 
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The true anomaly is given by 


02 



[1 — ecos (-E)]^ 


dM 


They can also be computed recursively using the Newcomb 
operators. 


5g2 

= M + 2e sin (M) H—— sin (2M) + e 


13 


+ e" 


103 . 11 . 

— sin (4M) - — sm (2M) 


12 

+ d (e®) . 


sin (3M) - - sin (M) 


(A3) 


Using the above expressions we can obtain the 
power series (in e) for (r\ 2 /cos (mij) — ^(j} 2 ) 
and (ri 2 /a)* cos (ttk)) — ^,(^ 2)1 ^-nd gathering 
cos [m(j) — (/r + n) UbI] terms [equations (IIOII and (1151) 1. we 
find the coefficients and 

Using computer algebra, we have computed „ and 
for jnj 8. Each coefficient is proportional, to leading 
order, to e'"', with higher order terms proportional to 

and so on. We have included terms up to order , so 
that every includes at least one higher order correction 

in eccentricity after its leading order. The leading terms for 
jnj ^ 4 are shown below. 

= -e [1 ± 2/i + 1] 

Ci '(^,±2 = + (5 ± 4/r)l + 4^^ ± 9/i + 4 ] 

C?X±z = + 6(2 ± + (12/i" ± 45/r + 38) I 

±8^® + 42/r^ ± 65/r + 27] 

^“±4 = 3^e'‘ [I* + (22 ± 8fi)L + (24ai^ ± 126/^ + 155) L 
+ (±32/r^ + 240^1^ ± 558^ + 390) I + ± 152fR 

+499fR ± 646/r + 256] 

(A4) 


^CB 

= 1 



^CB 

= [-1 ± 2ij] 



^CB 

= ie^ [L - (3 ± 4^)1 + ± 5)] 


^CB 

= + 3(3 ± 2/r)l' 

’-(12/.= 

± 33/. + 17) 1 


+2/r {±4fR + 15ii ± 13)] 



ryCB 

= [f-2{9±4^^)L 

+ (24/i^ 

± 102/r + 95) /■ 


-2 {±16fR + 96/r^ ± 165/r + 71) I 
+M (16/r® ± 120 m^ + 283/r ± 206)] 

(A5) 


Note that these coefficients are equivalent to the Hansen 
coefficients A“’^(e) (Murray & Dermott 1999), with 

and = x'^Je). (A6) 

The Hansen coefficients are defined according to 

exp [-1^02(0] = exp , (A7) 

J N 

which can be inverted, giving the following expression: 



cos {N[E-e sin (E)] - 
[l-ecos(E)]-(^+i) 


(A8) 
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